Hybrid Simulation Safety: Limbos and Zero Crossings

Physical systems can be naturally modeled by combining continuous and discrete models. Such hybrid models may simplify the modeling task of complex system, as well as increase simulation performance. Moreover, modern simulation engines can often efficiently generate simulation traces, but how do we know that the simulation results are correct? If we detect an error, is the error in the model or in the simulation itself? This paper discusses the problem of simulation safety, with the focus on hybrid modeling and simulation. In particular, two key aspects are studied: safe zero-crossing detection and deterministic hybrid event handling. The problems and solutions are discussed and partially implemented in Modelica and Ptolemy II

Higher harmonics of azimuthal anisotropy in relativistic heavy ion collisions in HYDJET++ model

The LHC data on azimuthal anisotropy harmonics from PbPb collisions at center-of-mass energy 2.76 TeV per nucleon pair are analyzed and interpreted in the framework of the HYDJET++ model. The cross-talk of elliptic $v_2$ and triangular $v_3$ flow in the model generates both even and odd harmonics of higher order. Comparison with the experimental data shows that this mechanism is able to reproduce the $p_{\rm T}$ and centrality dependencies of quadrangular flow $v_4$, and also the basic trends for pentagonal $v_5$ and hexagonal $v_6$ flows.Comment: 12 pages including 13 figures as EPS-files; prepared using LaTeX package for publication in the European Physical Journal

Modulational and Parametric Instabilities of the Discrete Nonlinear Schr\"odinger Equation

We examine the modulational and parametric instabilities arising in a non-autonomous, discrete nonlinear Schr{\"o}dinger equation setting. The principal motivation for our study stems from the dynamics of Bose-Einstein condensates trapped in a deep optical lattice. We find that under periodic variations of the heights of the interwell barriers (or equivalently of the scattering length), additionally to the modulational instability, a window of parametric instability becomes available to the system. We explore this instability through multiple-scale analysis and identify it numerically. Its principal dynamical characteristic is that, typically, it develops over much larger times than the modulational instability, a feature that is qualitatively justified by comparison of the corresponding instability growth rates

Solar wind interaction with comet 67P: impacts of corotating interaction regions

International audienceWe present observations from the Rosetta Plasma Consortium of the effects of stormy solar wind on comet 67P/Churyumov-Gerasimenko. Four corotating interaction regions (CIRs), where the first event has possibly merged with a coronal mass ejection, are traced from Earth via Mars (using Mars Express and Mars Atmosphere and Volatile EvolutioN mission) to comet 67P from October to December 2014. When the comet is 3.1–2.7 AU from the Sun and the neutral outgassing rate ∼1025–1026 s−1, the CIRs significantly influence the cometary plasma environment at altitudes down to 10–30 km. The ionospheric low-energy (∼5 eV) plasma density increases significantly in all events, by a factor of >2 in events 1 and 2 but less in events 3 and 4. The spacecraft potential drops below −20 V upon impact when the flux of electrons increases. The increased density is likely caused by compression of the plasma environment, increased particle impact ionization, and possibly charge exchange processes and acceleration of mass-loaded plasma back to the comet ionosphere. During all events, the fluxes of suprathermal (∼10–100 eV) electrons increase significantly, suggesting that the heating mechanism of these electrons is coupled to the solar wind energy input. At impact the magnetic field strength in the coma increases by a factor of 2–5 as more interplanetary magnetic field piles up around the comet. During two CIR impact events, we observe possible plasma boundaries forming, or moving past Rosetta, as the strong solar wind compresses the cometary plasma environment. We also discuss the possibility of seeing some signatures of the ionospheric response to tail disconnection events

Hamiltonian Hopf bifurcations in the discrete nonlinear Schr\"odinger trimer: oscillatory instabilities, quasiperiodic solutions and a 'new' type of self-trapping transition

Oscillatory instabilities in Hamiltonian anharmonic lattices are known to appear through Hamiltonian Hopf bifurcations of certain time-periodic solutions of multibreather type. Here, we analyze the basic mechanisms for this scenario by considering the simplest possible model system of this kind where they appear: the three-site discrete nonlinear Schr\"odinger model with periodic boundary conditions. The stationary solution having equal amplitude and opposite phases on two sites and zero amplitude on the third is known to be unstable for an interval of intermediate amplitudes. We numerically analyze the nature of the two bifurcations leading to this instability and find them to be of two different types. Close to the lower-amplitude threshold stable two-frequency quasiperiodic solutions exist surrounding the unstable stationary solution, and the dynamics remains trapped around the latter so that in particular the amplitude of the originally unexcited site remains small. By contrast, close to the higher-amplitude threshold all two-frequency quasiperiodic solutions are detached from the unstable stationary solution, and the resulting dynamics is of 'population-inversion' type involving also the originally unexcited site.Comment: 25 pages, 11 figures, to be published in J. Phys. A: Math. Gen. Revised and shortened version with few clarifying remarks adde

Developing an algorithm for pulse oximetry derived respiratory rate (RRoxi): a healthy volunteer study

Objective The presence of respiratory information within the pulse oximeter signal (PPG) is a well-documented phenomenon. However, extracting this information for the purpose of continuously monitoring respiratory rate requires: (1) the recognition of the multi-faceted manifestations of respiratory modulation components within the PPG and the complex interactions among them; (2) the implementation of appropriate advanced signal processing techniques to take full advantage of this information; and (3) the post-processing infrastructure to deliver a clinically useful reported respiratory rate to the end user. A holistic algorithmic approach to the problem is therefore required. We have developed the RROXI algorithm based on this principle and its performance on healthy subject trial data is described herein